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Abstract 

Models of modified gravity introduce extra degrees of freedom, which for 
consistency with the data, should be suppressed at observable scales. In the 
models that share properties of massive gravity such a suppression is due to 
nonlinear interactions: An isolated massive astrophysical object creates a halo 
of a nonzero curvature around it, shielding its vicinity from the influence of 
the extra degrees of freedom. We emphasize that the very same halo leads to a 
screening of the gravitational mass of the object, as seen by an observer beyond 
the halo. We discuss the case when the screening could be very significant 
and may rule out, or render the models observationally interesting. 
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1 Introduction and summary 



One of the most puzzling discoveries of our times is the fact that the present-day 
expansion of the Universe is accelerating [1]. Such an acceleration can be attributed 
to the existence of "dark energy"- a substance with a negative enough pressure - 
that is present in the Universe and, hence, in the rhs of the Einstein equation: 

G^, = StiGn{T^:'''' + T^t''" , (1) 

where G^jy stands for the Einstein tensor of the 4D space-time with metric gfj_^{x), 
and TJ^^^^^"^ and Tj^^"^^ energy (jg]2ote the stress-tensors for visible and dark matter, 
and dark energy, respectively. 

On the other hand, one can consider a different logical possibility: that the 
accelerated expansion is due to modified General Relativity (GR)|1. In this case, we 
would have the modified Einstein equations of the form: 

G^u - }C,A9, rUc) = SttG^T^T"'' , (2) 

where /C^,y((?,mc) denotes a tensor that could depend on a metric g, its derivatives, 
as well as on other fields not present in GR. Moreover, /C depends on a dimensionful 
constant rric ~ i^o ~ 10~^^ GeV, that sets the distance/time scale Tc = at 
which the Newtonian potential significantly deviates from the conventional one. For 
instance, in the DGP model /C^^^ is related to the extrinsic curvature tensor that 
gives rise to a self-accelerated solution [5l [6] (see, the comments on viability of this 
solution at the end of this section, and Ref. [7] for a recent review). 

Even though the difference between ([2]) and ([T]) might seem just conventional 
at a first sight, in reality, however, it could be observationally significant. For 
instance, it is possible to choose the rhs of ([1]) so that it gives rise to a background 
evolution obtained from ([2]), nevertheless, perturbations on these backgrounds would 
be different in ([2]) and ([1]), see, e.g., [8]- [IS]. 

In what follows we will focus on the issue of weather the two approaches, ([2]) 
and ([T]), could be differentiated by properties of a Schwarzschild-like solution for a 
static spherically symmetric source. 

A key feature of any theory of modified gravity of the form ([2]) is that, unlike GR, 
it allows for the possibility of having non-vanishing curvature outside of a source. 
This can be easily understood by taking the trace of ([2]) in a region outside of sources 
(where T^,^ = 0), that gives: 

-R = JC , (3) 

where JC = K,^ needs not be zero. This fact affects in particular the notion of mass 
and leads to the interesting phenomenon of screening. For example, when defining 

■^The latter approach is motivated by the "old cosmological constant problem" (for a review of 
which, see, e.g., [2], and in the context of modified gravity see, e.g., [3]). We will not be discussing 
this problem in the present work. 
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the Komar mass (which coincides with the ADM mass in stationary asymptotically 
flat space-time) one starts with an expression in terms of a volume integral of a 
projection of the Ricci tensor (see, e.g., [T7j) and then by use of Einstein's equations 
the Ricci tensor is replaced by T^,y — g^JT jl. In a modified gravity theory of the 
form ([2]) an extra term containing /C^j, and its trace is generated in the replacement 
of the Ricci tensor. Thus, the definition of the Komar mass contains an extra piece, 
referred to as a mass deficit below, which for a compact static source of spherical 
symmetry is given by 

AM oc Ml, j dr' (^/Coo - ^QoolC^ • (4) 

How significant is the mass deficit? The answer to this question would depend on 
a concrete model at hand. However, we would like to argue that in models which 
share properties of Lorentz invariant "massive gravity" , the mass deficit could be of 
the order of the mass itself. This has something to do with the fact that such models 
exhibit the so called strongly coupled behavior [T8],|19] (see also [201 [211 [221 [23] , and 
section 2 below for a summary), in spite of the fact that gravitational fields is weak 
everywhere [191 [211 [25l [26] . 




Figure 1 : Curvature extends outside a source to a distance r* . 

One way to interpret this property for a static spherically symmetric source is 
to observe [27] that the source gives rise to the curvature that extends up to a 
macroscopic distance r,,, in the way depicted in figure [H This curvature contributes 
to the integral 

In a concrete case of the non-perturbative Schwarzschild solution of the DGP 
model [27], we will show that this integral is saturated around a distance (~ 
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{M/mlMpiy^^ - the topic of the next section ) where the value of /C is typically of 
order m^. Then, the result of (jlj) is 

AM ~ M , 

i.e., the contribution to the mass from the modification of gravity term is of the 
same order as the mass itself! 

In the following sections we will make concrete the statements pointed out above 
in the DGP model of modified gravity. They include the perturbative arguments 
leading to the derivation of the scale and details on the exact solutions available 
for Schwarzschild-like sources and domain walls. While the former are based on 
an ansatz, legitimate concerns about the bulk boundary conditions of which were 
raised in [28], nevertheless, the fact that the ansatz recovers very precisely the 4D 
GR results at short distances and smoothly interpolates to the 5D regime (that 
no other solution is known to do), suggest that it may be capturing right physics. 
Moreover, the above properties were subsequently found to be true for the case of a 
domain wall for which an exact solution was obtained [29] . 

Although the existence of the mass deficit (jl]) could be interesting observation- 
ally, it may lead in certain cases to problems with the theory. Indeed, because of 
the terms (jl]) in the expression for a gravitational mass the proof of the positive 
energy theorem [301 EI] is not directly applicable. Hence, in general, there could 
exist negative "mass" solutions [SI EI] even for matter stress-tensors that satisfy 
conventional positive energy conditions. One example of this is the self-accelerated 
solution [3 dH]. Small perturbations about this solution in a linearized theory and 
with non-conformal sources exhibit ghost-like states [21], [231 Ell [331 El], however, 
there exist serious arguments that the perturbative results cannot be trusted in the 
full non-linear theory [53] (see also [3S])- Nevertheless, some semi-exact [27] and 
exact [221 [331 [3Z] non-perturbative solutions on the self-accelerated branch exhibit 
"negative mass" . This suggest that the self-accelerated branch should be unstable, 
however, it is not clear what is the time of its instability. An explicit calculation 
on decay of the selfaccelerated branch into the conventional one shows that such 
a decay does not take place, at lest in a quasi-classical approximation [3B]- This 
question is still open and we will not be discussing it further in the present work. 

If the mass screening is substantial, then, at scales beyond r^, gravity would be 
modified significantly. However, the value of the scale for the entire observable 
Universe is Hq^ ~ 10^^ cm. Therefore, on average, the beyond-r^,-effects will be 
hard to detect. There may be exceptions for isolated clusters of galaxies separations 
between which are greater than their own scales, and any other r^, scales in their 
vicinitj0. For precise calculations of the beyond r^, physics, however, new averaging 
technique and non-perturbative calculational methods would be needed. 

^We thank Lam Hui and Roman Scoccimarro for discussions on these issues. 
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2 The r*— scale (Vainshtein scale) 



In what follows we will concentrate on the concrete example provided by the DGP 
model [1] in which all interactions except gravity are confined to a 4D brane embed- 
ded in an infinite volume 5D empty space where gravity propagates. 

The modification of gravity in this model is given in terms of the extrinsic cur- 
vature K^i, of the brane and reads: 

= mc{K^^ - g^uK) . (5) 

The 5D space has coordinates {x^,y) with the brane at the surface y = and the 
4D metric in ([5]) is g^u{,x^, 2/ = 0). 

The macroscopic distance in the DGP model can be derived by considering the 
linearized analysis of the theory j^. In particular we will focus on the one-graviton 
exchange amplitude between two sources whose gauge independent expression reads 
as follows: 

^l-graviton(P, V) = o exp{-p\y\) , (6) 

where 

Ti% = SttG^It'^^ - ^-T-T^ , (7) 

and is the square of the Euclidean brane 4-momentum. The pole at = has 
zero residue, and the second pole in ([H]) is on a non-physical Riemann sheet. The 
former implies the absence of a massless graviton in the exchange while the latter 
describes the propagation of a metastable state with lifetime ~ m~^, which decays 
into a continuum of KK modes. 

The striking feature of ([7]) is that in the rric — ^ limit (in which the modifica- 
tion of gravity should disappear) the numerator does not reduce to the analogous 
expression in GR: 

SvrG^^T^^ - ^T-t) . (8) 

This fact, which could be used to exclude (I7j) by observations, is known as the van 
Dam-Veltman-Zakharov discontinuity (vDVZ) [39]. The difference between and 
([8]) is due to the fact that a 5D graviton (or a massive graviton for that matter) 
propagates 5 on-shell degrees of freedom (helicity-2, helicity-1, and helicity-0), while 
the GR graviton propagates only 2 on-shell degrees of freedom (helicity-2 state). 
And while the helicity-1 state of the 5D graviton does not contribute to at the 
linearized level because of the contraction with conserved sources, the helicity-0 
state couples to the trace of the energy- momentum tensor and gives a non- vanishing 
contribution (when T ^ 0). 
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It has been argued in [181 [IH] that the continuity in the — hmit would be 
restored if nonhnear effects were taken into account. The relevance of these terms 
can be understood in the following terms: the longitudinal part of the graviton 
propagator in DGP contains terms proportional to PuPiy/nicp. This term does not 
contribute to the amplitude because of conservation of the stress-tensor, but it 
does contribute already in the first nonlinear correction (since the stress-tensor is 
only covariantly conserved in the non-linear theory). And due to the singular behav- 
ior of these terms in the —>■ limit, perturbation theory breaks down prematurely. 
However, this breakdown is an artifact of an ill-defined perturbative expansion - the 
known exact solutions of the model have no trace of breaking [19]. The perturbative 
expansion in powers of gets "contaminated" by another dimensionful parameter 
l/rric, and this leads to its breakdown. As possible ways forward one could either 
adopt a different type of expansion, e^., an expansion in the small parameter rric 
[T9| [21] , or try to find exact solutions U- Both of these programs have been carried 
out to a certain extent and we will review in the next section the salient features of 
the latter. 

As presented in the DGP model has one adjustable parameter, namely rric 
which determines a scale that separates two different regimes of the theory. For 
distances much smaller than one would expect the solutions to be well approxi- 
mated by GR and the modifications to appear at larger distances. This is indeed the 
case for distributions of matter and radiation which are homogeneous and isotropic 
at scales > Tc [5], [6], [19] . However, more compact sources exhibit different proper- 
ties. For example, a compact static source of the mass M and radius tq, such that 
i^M < <^ Tc {tm = 2GnM is the Schwarzschild radius) a new scale, combination 
of Tc and rM, emerges (the so-called Vainshtein scal^) [T9] : 

r^^ivMrl)'/'. (9) 

Below this scale the predictions of the theory are in good agreement with the GR 
results and above it they deviate considerably. These type of sources will be dis- 
cussed in more detail in the next section together with other ones with higher simpler 
symmetry: domain walls. 

3 Concrete examples 

In this section we will focus on the mass screening phenomenon describing how 
it arises in the cases of the Schwarzschild-like non-perturbative solution (NFS) of 
[271 US] and the exact domain wall (DW) solutions of [29] . 

^It is also possible to modify the theory at the linearized level so that the conventional pertur- 
bative expansion is well-behaved [45l l46] . [47] . |48j . 

similar, but not exactly the same scale was discovered by Vainshtein in massive gravity |18j . 
hence the name. 
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3.1 Schwarzschild solution 



The NPS solution studied in [271 US] is found by considering a static metric with 
spherical symmetry on the brane and with Z2 symmetric line element: 

ds^ = -Q-^dt^ + e^rfr^ + r'^d^f + 27 drdy + e^dy"^ , (10) 

where A, 7, a are functions of the radial coordinate on the brane r and the transverse 
direction coordinate y. The Z2 symmetry across the brane {y = 0) implies that 7 is 
an odd function of y while the rest are even. The choice —gu = 1/grr represents an 
ansatz, but notice that we have kept the off-diagonal term Qry = 7. 

The brane is chosen to be straight in the above coordinate system but one could 
transform ffTOj) to another one in which the metric is diagonal ds"^ = —A{r, z)dt^ + 
5(r, z)dp^+C{r^ z)dVt^+dz^. Here, 1/5 in general, and the ansatz is reflected in 
the fact that in this system our brane will be bent or, in other words, in a particular 
nontrivial choice of the position of the brane z{r). 

This ansatz allows us to close the system of equations on the brane, and leads to 
the following solutions for the gravitational potential 0, 



20 , (11) 
dr r^U{r) . (12) 



1^ N , ^ 2 

(h(r) = ? ' -^^ - 

The function f/(r) in (fT2|) is given implicitly by the solutions of the following two 
equations (giving rise to a conventional and self-accelerated branch respectively): 

t/2(3 + 3t/ + V3/)2v^(-5 - 3t/ + /) ' 

(kr)' - (-5-3^ + /)(-3-3^-V3/)^-/^ 
- (f/ + 2)2(1 + 3t/ + /) ' 

where / = y/TTWTSlP and k is an integration constant. 

The off-diagonal and yy metric components are determined from 



(r7e" 



(r^7e 



IT 

-\ 2 



(15) 



+ e-V, (16) 



and the profile of the warp factors ( "y-derivatives" ) can be computed on the brane 
as well. 

There are two integration constants, k and the one produced in the integration 
( IT2l) . which are determined by imposing appropriate boundary conditions near the 
source (r <^ and at large distances. For the first condition we impose the 4D 
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behavior of the potential near the source: = G^M/r, while for the second one 
we require that the coefficient of the possible 1/r term be zero, i.e., to be left with 
5D behavior at large distances, namely, A ~ ^'m/^'^ in the conventional branch and 
A ~ m^r^ + f^i/r"^ in the self-accelerated branch. 



3.1.1 Conventional branch 

The conventional branch is obtained from the solution of (fT3|) . As shown in [2j 
the boundary conditions discussed above determine the asymptotic behavior of the 
solution. At short distances, r <C (f/ — > +00), we get 

(17) 



r 2 V r 

where ai ~ 0.84 and the coefficient of the 1/r term was chosen to be —G^M by 
fixing the constant of integration in (1121) . 

The other integration constant ki is chosen such that at large distances, r ^ 
{U —>■ O"*"), we obtain an expansion with no 1/r term: 

~2 

^ - + . . . , (18) 



2r2 

which fixes the value of ki in terms of r* ((r*fci)^ ^ 0.21) and also implies 

fl,^ ^ 0.56 TMr, . (19) 

The relation (fT9|) should be contrasted with the naive expectation from a lin- 
earized analysis: in the 5D regime one would have expected to have ~ fM^c, 
however, we get a much smaller value, reduced by a factor r^/vc = ij-j^i /r^^/^ . 

Therefore, as we see, a short distance observer at tm r <^ would measure 
the gravitational mass M with a small corrections to Newton's potential, while the 
large distance observer at r 3> r* would measure an effective gravitational mass 
~ M{rM /"TcY^^ [27]. The latter includes the effects of the 4D curvature which is 
significant up to r^, as depicted in figure [2l The 5D mass is partially screened at 
large distances. 

Another point worth emphasizing is that a perturbative expansion suggests that 
for r* -C r -C rc the metric should have an approximately four-dimensional, 1/r, 
scalar-tensor-gravity type form [ll[T9]. However, the NPS above exhibits a different 
behavior: beyond r,, the metric turns into the one produced by a five dimensional 
source. We interpret this as a complete screening of the 4D mass of the source by 
the halo of non-zero curvature. 

The screening of the 4D mass can be made explicit by taking into account the 
expression for the 4D Komar mass as a function of r. In the ansatz used here gives 
this gives an effective mass 

Meff ^ Mli{-r(t) + T{r(t))r) . (20) 
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Figure 2: Conventional branch 



The first term the rhs of fl2U]) is a smooth decreasing function of r and gives a 
contribution that is ~ M (the original mass of the source) up to r^, and rapidly falls 
off like 1/r beyond that point. The second term, on the other hand, depends on the 
gradient of and is peaked around r^. The combined effect is seen in figure [S] the 
effective mass increases from its 4D value M near the source up to r ~ and then 
falls to zero abruptly. 

1.5 . ■ . ■ — — 1 




0.3 ' ■ ' ■ — ' ■ — 

0.150.2 0.3 0.5 0.7 1 1.5 2 



Figure 3: Log-Log plot of M^ff/M vs. r/r^,. The oblique dashed line shows a 1/r fall-off. 

The short distance mass increase can also be deduced form the approximate form 
of the potential (fT7|) since the second term provides an additional attraction toward 
the source. 
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3.1.2 Self-accelerated branch 



The solution on the self-accelerated branch is obtained from ( fl^ . The relation 
between ki and is obtained, as in the conventional case, by imposing boundary 
conditions together with the following asymptotic behavior. At large distances, 
r r^: {U we derive 

where, 

rl^ ^ 0.45 TMn , (22) 
, i.e., 5D mass screening. At short distances, r ^ r* (f/ — > — oo), we get 

GnM 1 2 2 /r A 2(^-1) 



where Q!2 = — ai ~ —0.84 is, in absolute value, the same constant appearing in the 
conventional branch short distance expansion (fT7|) . Note, however, that the sign of 
the correction to the 4D behavior is opposite in the two branches. 




Figure 4: Self-accelerated branch 

At intermediate distances, r^. ^ r ^ r^, the potential contains a 5D gravitational 
term that is repulsive, f\Jr'^. This looks like a 5D negative mass. However, this 
is not an asymptotic value of the mass since one can only cover the solution in 
the above coordinate system till r ~ Tc where the dS like horizon is encountered. 
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Moreover, in the intermediate regime r^, -C r <^ r^, the de Sitter term m^r^ in the 
potential always dominates over the f\j/r'^ term suggesting that the effects due to the 
Schwarzschild source are strongly suppressed. The picture that explains screening 
in this branch is depicted in figure HI 

3.2 Domain Walls 

The second example that illustrates the screening phenomenon in the DGP model 
is that of a brane DW source. The study of this type of sources was done in [29] on 
which we base the following discussion. 

The source considered is a Nambu-Goto DW in 4D localized at 2; = with 
stress-tensor 

T^, = a (5(^)diag(l, -1,-1,0) , (24) 

where a is the tension of the wall and z denotes the coordinate transverse to its 
world- volume spanned by {t,x,y). 

The domain wall solution for such a source in GR [l9l [50] displays 3D de Sitter 
expansion in its world- volume at a rate H = 27rGNO'. 

In DGP, however, the situation is different. For tensions smaller than a critical 
value, (Tc = = ■mcMpi/2, the wall has no gravitational effects. One way to 
understand this is by noticing that for these sources, the modification of gravity 
term ([5]) precisely compensates the energy momentum tensor T^j,. Therefore, the 
tension of the wall, as seen from the point of view of a 4D observer, is screened 
entirely by gravitational effects encoded in the extrinsic curvature. Not surprisingly, 
the domain wall world- volume remains fiat, and so does the metric on the brane. 

Furthermore, this screening takes place inside the core of the wall. Hence, the 
analog notion of the scale for a domain wall (understood as where the self-shielding 
takes place) coincides with its thickness, 

r^^^^ = d 

' ^ '^core • 

This is to be compared to the Schwarschild-like case, where the shielding also oc- 
curs, and extends outside the source. The net result is the screening of the 4D 
tension/mass in both cases. 

For supercritical tension branes the extrinsic curvature can no longer balance 
the energy momentum tensor and the brane inflates. The transverse direction to 
the wall is compactified to a size d and a zero mode graviton appears. Since this 
phenomenon also takes place in GR, it provides a means to contrast the 5D effects 
in DGP with those of a supercritical DW in 5D, i.e., by comparing the world- volume 
inflation rate in both cases. 

The exact solutions of [22] give a suppression factor for the inflation rate of 
the supercritical DW in DGP for both branches of solutions. In the conventional 
branch it is d/2rc while in the self-accelerated branch it is given by —d/2rc. As 
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argued above, we should identify the Vainshtein scale for these sources d and 
therefore we find the same parametric screening of the 5D tension as for the mass 
in the Schwarzschild-like solutions of the previous example. 



4 Discussions 

The above results may also be applicable to other models where the results of 4D 
gravity are recovered through strongly coupled behavior. The minimal model of 
brane induced gravity in greater that five dimensions [51] contains ghosts [52] . 
However, its extensions are ghost free [53], a small subset of which has a strong 
coupling regime [MIIMI (see, also [55]). A recent model of cascading brane induced 
gravity is also ghost free [56]. It would be interesting to understand the issue of 
presence/absence of the mass screening in these models (see, also, [57]). 

Unfortunately, at present there is no consistent 4D theory of Lorentz invariant 
massive gravity, as it suffers from nonlinear instabilities [58l [59| [60l [6T] . However, 
the mass screening effect described above is based on rather universal principles, 
and it is reasonable to expect that the phenomenon will also persists if a consistent 
model is found. 

It would be interesting to understand the issue of the mass screening in the 
/(i?)-type models of models of modified gravity, see, e.g., [62], as well as in Lorentz 
violating models [Ml Ell ESI ESI EI] . Some of these models [Ml ES, EEl E7] avoid the 
strong coupling behavior, in which case the mass screening is not expected to be 
very significant. 
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